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Abstract 



We present a gauged Lifshitz Lagrangian including second and forth order spatial derivatives 
of the scalar field and a Chern-Simons term, and study non-trivial solutions of the classical 
equations of motion. While the coefficient /3 of the forth order term should be positive in order 
to guarantee positivity of the energy, the coefficient a of the quadratic one need not be. We 
(~| investigate the parameter domains finding significant differences in the field behaviors. Apart 

from the usual vortex field behavior of the ordinary relativistic Chern-Simons-Higgs model, we 
find in certain parameter domains oscillatory solutions reminiscent of the modulated phases of 
^ Lifshitz systems. 
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1 Introduction 



In their well-honored proposal to describe dual strings [T] , Nielsen and Olesen stressed the connection 
between the Abelian Higgs model and the Ginzburg-Landau theory of superconductivity, relating 
the free energy in the latter with the action for static configurations of the former. In this way, the 
vortex filaments of type-II superconductors were identified with string-like classical solutions in a 
gauge theory with spontaneous symmetry breaking. 

More than 30 years ago Ginzburg proposed [2] a generalization of the Ginzburg-Landau func- 
tional by including higher derivative terms, this implying an anisotropic coordinate scaling, in order 
to describe superdiamagnets - a class of materials with strong diamagnetism but differing from 
conventional superconductors. Such generalization was then used to analyze [3] the properties of 
superconductors near a tricritical Lifshitz point, a point in the phase diagram at which a disordered 
phase, a spatially homogeneous ordered phase and a spatially modulated ordered phase meet. 

The study of Lifshitz critical points has recently attracted much attention, not only in con- 
nection with condensed matter systems (see |3] and references therein) but also in the analysis of 
gravitational models in which anisotropic scaling leads to improved short-distance behavior (see f5] 
and references therein). A link between these two issues was established in [6j within the framework 
of the gauge/gravity correspondence by searching gravity duals of nonrelativistic quantum field 
theories with anisotropic scaling, dubbed in [8] as "Lifshitz field theories". 

The question that we address in this work is whether one can find Nielsen-Olesen like solutions 
when anisotropic scaling is introduced in the Abelian Higgs model through the addition of higher 
order spatial derivatives. As a laboratory we consider a 2 + 1 dimensional model with a complex 
Higgs scalar coupled to a U{1) gauge field with a Chern-Simons (CS) action [9]. The topological 
character of the CS term avoids the possibility of including higher order derivatives for the gauge 
field action (as it would be case for the Maxwell action). 

When higher order derivative terms in the scalar Lagrangian are absent, the Chern-Simons-Higgs 
model has vortex-like finite energy solutions carrying both quantized magnetic flux $ and non trivial 
electric charge Q = — k<I> with k, the CS coefficient |10|-|11). Moreover, for an appropriate sixth- 
order symmetry breaking potential, ffist order BPS equations |12)-|14j exist, which can be easily 
found by analyzing the supersymmetric extension of the model |15| . Our goal will be to determine 
whether this kind of solutions also exists in a "Lifshitz Abelian Higgs" model and, in the affirmative, 
how they depend on the parameters associated to the Lagrangian scaling anisotropy. 

The plan of the paper is the following: we introduce in section 2 a {2 + l)-dimensional Lifshitz- 
Higgs model with gauge field dynamics governed by a Chern-Simons term. In order to solve the 
classical equations of motion we make the same ansatz leading to vortex solutions in the ordinbary 
(relativistic) case. Then, in section 3 we analyze the asymptotic behavior of the gauge and scalar 
fields resulting from the equations of motion, showing the existence of four regions according to 
the values of the parameters of the model. We discuss in section 4 the properties of the solutions 
obtained numerically i giving a summary of results and a discussion on possible extensions of our 
work in section 5. We briefly describe in an appendix the linearized approximation we employed to 
determine the asymptotic behavior of the solutions in different parameter regions. 

2 The Lagrangian 

We consider a 2 -|- 1 dimensional model with Chern-Simons-Higgs Lagrangian 

L = 7 \Do[A]cl)\^ - a \Di[A]<p\^ - /3 \Di[A]Di[A]^\^ + + |e/^--^^9,A„ (1) 
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with fi = 0,1,2 and i = 1,2. The metric signature is (1,— 1,— 1). We consider space and time 
coordinate units so that 

[x? = [t] . (2) 

Accordingly, 7, /3 and n are dimensionless and a has length dimensions [a] = —2. Concerning the 
dimensions of the complex scalar (j) and U{1) gauge field one has [(/>] =0, [Ai] = — 1, [Aq] = —2. 

The Lagrangian ([T]) is a generalization of the one considered in [12J-[13] incorporating higher 
(forth) order covariant derivative terms for the scalar fields. For vanishing potential and at the 
"Lifshitz point" a = 0, the Lagrangian is invariant under anisotropic scaling with "dynamical critical 
exponent" z = 2 

Xx , \h . (3) 

Note that the choice of a Chern-Simons term ensures that scale invariance is preserved even in the 
presence of gauge fields (as opposed to what would happen with a standard Maxwell term). 
The covariant derivative acts on the scalar field (j) according to 

D^[A](t>={d^ + ieA^)<l, (4) 

with [e] = 0. The potential V[(j)\ is to be specified below. 

Given the Lagrangian ([T| one gets Gauss's law by differentiating with respect to Aq, 

ne'^'^d.Aj = f (5) 

where 

jo = ie^{4>*Do(t> - (pDocI)*) = -2e Vo|<^|' (6) 

Defining 

B = -e'^diAj (7) 

one then has, using eq.([5]), 

K B 
^'=2^W 

Inserting this result in eq.([6]) one gets 

JO = (9) 
so that the usual Chern-Simons-Higgs model relation between charge Q and magnetic flux $ holds 

Q = j d^xf = j d^xB = -K^ (10) 

The energy density £ associated to Lagrangian M is 



£ = a mA]ct>\^ + /? I A[^] A[^]<Al' + ^^ + ■ (11) 

A lower bound for the energy requires (3 to be positive while a can have any sign. 

As stated before, in the /3 = 0, 7 = 1 relativistic case and for a sixth order symmetry breaking 
potential this theory is known to have, at the classical level, self-dual vortex solutions both in the 
Abelian case |12)-|13) and in its non-Abelian extension |14) . 

In order to solve the Euler-Lagrange equations deriving from Lagrangian ([T| we consider the 
static axially symmetric ansatz 



4> = f{r) exp{—irup) (12) 

Mr) , , 

A^ = ^ 13 

r 

^0 = Ao{r) (14) 
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with n G Z. Given this ansatz the magnetic and electric fields read 

r dr ^ dr 



, 1 dA(r) ^, , dAn , . 

B{r) = , E{r) = — (15) 



The equations of motion take the form 
K dA{r 



+ 27eMo(r)r(r) = (16) 



dr 

o2 . -jX o„2 



^dvlo _^ 4e2/3 /n ^ / _^ 1 d ^ j ^ + A] f - {) (17) 

dr r V e / \ dr"^ r dr r^Ve // rVe / 

\ dr"^ r dr V e ^ / \ '^'''^ dr V e / 

fd^f ^idf 1 2 \ 2 2 lay 



\ dr^ r dr r"^ J 2 9/ 

The potential is in general chosen so as to allow for spontaneous symmetry breaking. In the 
relativistic 2 + 1 dimensional case the most general renormalizable self-interacting scalar potential 
is sixth order and in fact to find first order BPS equations it should be of this order and take the 
form p]-[T3] 

y = P,nf-v'? (19) 



with V the Higgs field vev and the coupling constant r has length dimensions [r] = —2. In the 
relativistic model first order self dual equations exist at a certain value r = tbps which would 
correspond in the present Lifshitz case to tbps = From here on, and in order to compare 

the Lifshitz model results with those arising in the relativistic case, we shall take V as given in (19) 
and r = tbps- 



3 Asymptotic behavior 

We start by discussing the conditions that we shall impose at the origin and at the boundary. We 
choose as conditions at the origin those leading to regular solutions in the relativistic case (see for 
example [12]): 

f{r) = /orl"l 
Ao{r) = ao + cor^l"! r 

A{r) = dor2|"l+2 (20) 

Note that a constant term oq in the ^o('') expansion is included in order to achieve consistency of 
eq.(16) at the origin. Coefficients gq and do are related according to 

Concerning large r, we write 



/(r) PS v + h{r) (22) 

71/ 

A(r) ha(r) r ^ oo (23) 

e 

Aoir) ^ ao(r) (24) 
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with h{r), a(r) and ao{r) small fluctuations. We then linearize the equations of motion which reduce 
to 

- /3 V^V^ h{r) + aVl h{r) - ah{r) = (25) 
1 da{r) 



r dr 

dao{r) afi 



dr r 

where 



+ 7/iao(r) = (26) 
a(r) = (27) 



dr"^ r dr ' 



Eqs.(26 )-(27) can be written as two decoupled second order equations 



dr"^ r dr 

First we deal with the scalar field behavior. After writing 



— ^ H -, a7/i ao = (30) 

ar^ r dr 

— ajfi'^a = 0. (31) 



d'^a 1 da ^ 



^(^) = 6xp(gr), (32) 



with a constant, the solutions are determined from the equation 



4 = ^(a±Va2_4^^y (33) 



The asymptotic behavior of the scalar field is then given by 



ir 

From the results one can see that there is a critical value for /? 

4C7 



/(r) w + — exp(-g±r) (34) 



/3crit = ^ (35) 



above which become imaginary. 

In the region q > and /3 < /3crit the solutions for q± are real. In particular, for fia <C 

'Z--- (36) 

Note that g_ coincides with the standard relativistic case solution where it plays the role of the 
Higgs field mass [I2j. 

Concerning the region a > and /3 > fi^rW one has 



d = 5^ (37) 
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which gives a complex solution. This region corresponds to underdamped oscillations of the Higgs 
field. We can write this as 



ql = ^e^^^ (38) 



where 



The solution is therefore 



where 



a 



2 



tan(x) = ^ . (39) 

a 

cos{k r + 6) (40) 



a 

exp (— Ar) 



A = ^^1 cos(x/2)| , k = ^^\ sin(x/2)| (41) 

where 5 is a constant phase. 

We now consider the case of q < 0. In this case for f3 < /3crit we have that 

Q± = ^ (-I«l ± \/«'-4a/3) (42) 

which is always negative leading to oscillatory solutions with wavenumbers \q±\. 
Finally let us consider the /3 > /3crit region where the solutions become 

ql = ^e^^^ (43) 

leading for the scalar field behavior to a situation similar to the case of a > with /3 > /3crit- 

Let us now study the asymptotic behavior of the gauge fields. For a > the consistent asymp- 
totic behavior is 

ao(r) w -^exp(— Air) 

a(r) aoo\/r exp{—kr) (44) 

Notice that in this region the asymptotic field behavior ensures finite energy and quantized magnetic 
flux as in the relativistic case 

$ = — n , n e Z (45) 



In the a = case linearization leading to eqs.(30)-(31 ) is no longer valid. Instead, writing 
a = y/rg{r) and using the gauge field equations of motion one gets a second order nonlinear equation 
for g compatible with bounded solutions at infinity. As will be discussed in next section, we do find 
a bounded numerical solution for a = 0. 

Concerning the a < region, one has 

ao{r) ^^^^ sin(A;r + (^) 
a{r) ^ aoo\/r cos(A;r + (^) (46) 

with 

k = v^|a|7/i , Ooo = - W 7^ aooo (47) 

V a 



The oscillatory behavior of configurations satisfying (46) will require the introduction of appropriate 
boundary conditions at a finite radius R. 
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4 Solutions 



We shall present in this section numerical solutions of eqs.(18| satisfying the asymptotic condition 
discussed above. For definiteness we take n = 1 and we shall fix 7 = 1 (Since we are considering 
static solutions, changing gamma amounts to a redefinition of the scalar field coupling with Aq). In 
order to ensure positivity of the energy we shall take /3 > 0. We shall separately consider a > and 
a < regions. Following the discussion in the previous section, we shall distinguish regions with 
/3 ^ Pcrit- The numerical procedure is based on a forth-order finite differences method applied in 
the interval (e, R) with e close to the origin and R large, in combination with the behavior of fields 
close to the origin given by eq. ( 20 ) . 



4.1 The a > region 

We start by studying the the a > 0, /3 < /3crit region. We give the results of our numerical calculation 
for E and B in figures 1 and the scalar field in figure 2. 




Figure 1: The electric (solid line) and magnetic (dashed line) fields in the region a > 0, /3 < f3crit, 
with a = 1, /3crit = 0.0625 and /3 = 0.04. As in the relativistic Chern-Simons-Higgs model, the 
magnetic and electric fields form a ring surrounding the vortex core. 



One can see that the profile of the fields in this region exhibit slight deviations to the relativistic 
case, originated by the fourth order derivative terms. It should be noted that as /3 grows we found 
numerically that the maximum magnitude of the electric and magnetic fields decrease. Concerning 
the Higgs field, it reaches its vacuum value exponentially according to eq. (34), as can be seen in 
figure 2 with a similar profile as that corresponding to the relativistic case as shown in figure 2. 
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Figure 2: The Higgs field profile in the region corresponding to a > 0, /3 < /3crit- (a = l,/3crit = 
0.0625 and /? = 0.04) 



Let us now consider (3 > /3crit range where the roots q± are complex (37), this giving rise to 
underdamped oscillations in the Higgs the profile, as shown in figure (3). 




2 4 

Figure 3: The Higgs field profile in the region a > 0, /3 > /3crit- We have chosen a = 1, /3crit = 0.0625 
and /3 = 0.2. The inset shows a zoom of the region where / overshoots its vev and comes back to 
it, as is characteristic of an underdamped behavior. 



For a given value of a the magnetic and electric field solutions for (3 > ficrit are qualitatively 
the same as those shown in figure 1 for /3 < Pcrit- 

We then conclude that in the a > region the electric and magnetic field behavior is very 
similar to the ordinary relativistic CS-Higgs model. Concerning the scalar field, as one crosses 
from /3 < /3crit to /3 > /3crit, it changes from the usual to an underdamped approach to its vacuum 
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expectation value. 

We have studied the /3-dependence of the energy in this region finding a hnear behavior for small 
f3. As an example, we show in figure 4 a numerical calculation of the energy £ as a function of f3 
for a = 1, /Scrit = 0.0625. We find that £ behaves approximately as £ ^ £q + 0.25/3. 
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Figure 4: The energy as a function of /? for a = l,/3crit = 0.0625 



We end this subsection by discussing the a = case for which, for vanishing potential, the 
Lagrangian is invariant under anisotropic scaling with "dynamical critical exponent" z = 2. In this 
case /3crit = so that for any (3 > the Higgs field shows an underdamped behavior. We have 
numerically confirmed this result and also found bounded solutions for the gauge fields. The field 
profiles are qualitatively similar to those found for a > 0, /3 > PcrU- 



The a < region 



One expects in this region a clearly different behavior compared to the relativistic CS-Higgs system 
since the negative sign of a in the \Di(j)\'^ energy term implies not only a change of sign in the |V(/>p 
term but also in the gauge field "mass term" that now has the "wrong" sign. 

We start by studying the /3 > /3crit region where the fields asymptotic behavior is given by eqs. 



(43)-(47). This behavior leads to an oscillatory energy density (and consequently to an in general 



unbounded energy). For example, the third term in expression (11 ) for the energy density takes the 
asymptotic form 

We show in figure 5 the electric and magnetic fields in the a < 0, P > /3crit region. Their 



profiles show the asymptotic oscillatory damped behavior consistent with eq. ( 46 ) . The behavior of 



the scalar field is presented in figure 5. A zoom outside the vortex core shows damped oscillations 



consistent with equations (34)- (43) 
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Figure 5: The electric (solid line) and magnetic (dashed line) fields in the region a < 0, /3 > f3crit, 
with a = -0.2, /3crit = 0.0025 and /? = 0.25. 




2 4 6 8 10 

Figure 6: The Higgs field profile in the region a < 0, (3 > (icriti with a = —0.2, /3crit = 0.0025 and 
/3 = 0.25. 



In the ;S < /3crit region, the roots we found in section 3, eq.(42 ), lead to pure oscillatory solutions 
with no damping. The assumption of h in eq.(22) being asymptotically a small perturbation to 
the scalar vacuum expectation value v is then not self-consistent. We have not been able to find 
stable solutions of our 2 + 1 model with the ansatz (12)-(13). We indeed know that in the absence 
of dynamical gauge fields this range of parameters corresponds to the modulated ordered Lifshitz 
phase associated to spontaneous breaking of translations p!9]. We then conclude that in this region 
a more detailed numerical study allowing the implementation of more general ansatze would be 
necessary. 
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5 Summary and Discussion 



We have proposed a gauged Lifshitz Lagrangian with higher (forth) order spatial derivatives of the 
scalar field and a CS term and studied numerically non-trivial solutions of the classical equations 
of motion. Notice that contrary to previous analysis of Lifshitz theories with CS term [21j with 
z = 2 we considered higher derivatives for the scalar field rather than for the gauge fields. As a 
consequence, the classical solutions of our model have a different character of the ones resulting 
from such model |22j . 

Coming back to the model we analyzed, let us recall that /3, the coefficient of the forth order 
derivatives term, was taken positive in order to ensure positivity of the energy. In contrast, the 
a coefficient multiplying the ordinary second order derivative term could take both positive and 
negative values being a = the Lifshitz point at which the model exhibits z = 2 anisotropic scaling 
in the absence of a potential term. 

In order to solve the equations of motion we have made the static axially symmetric ansatz that 
leads to vortex solutions in the relativistic case. For a > we have found solutions with magnetic 
and electric fields qualitatively similar to those of the ordinary relativistic model. The magnetic 
flux is quantized and the usual relation between electric charge and magnetic field in CS systems 
holds. The difference with the standard relativistic case manifests more pronouncedly in the Higgs 
field behavior which for /? > f^crit approaches its vacuum expectation value with underdamped 
oscillations. The critical value is given by formula (35), ficrit = o?l^cP'^ showing a dependence 
on the coefficient of the quadratic derivative coefficient and on the parameters of the model (the 
value V of the Higgs field at the minimum, the gauge coupling e, the CS coefficient k and the Higgs 
field self- interaction coupling constant r). For a = the numerical solutions that we found are 
qualitatively similar to those found for a > 0,/3 > /3crit- 

The situation for the a < region radically changes basically because of the change in sign of 
the gauge field mass term. The ansatz led to pure oscillatory solutions for the gauge fields with no 
damping. Concerning the scalar field one can again distinguish two situations depending on wether 
fi is larger or smaller than /3crit- In the former case we were able to find solutions exhibiting electric 
and magnetic field profiles with an asymptotic oscillatory behavior while the Higgs field profile 
shows damped oscillations. This behavior leads in general to an oscillatory energy density and an 
unbounded energy. In the /3 < (icrit region, the proposed ansatz led to pure oscillatory solutions 
with no damping. 

We think that in the region a < other terms in the Lagrangians, as those considered by 
Ginzburg for the free energy of superdiamagnets and superconductors [2j might become relevant. 
Also, more general ansatze, not purely relying in cylindrical symmetry should be considered in 
order to incorporate the possibility of asymptotic breaking of translational symmetry which is 
characteristic of modulated Lifshitz phases. We hope to come back to this problem in a future 
work. 



Appendix: The Frobenius Method 

In this section we wish to apply Frobenius's method to the linearized Higgs field equation of motion 



in order to determine its behaviour close to r = 0, where / is assumed to be small (see eq. (20)) and 
the equation has a regular singular point. Following ref. [2^ we recast the equation of motion (17) 
for the Higgs field close to r = in simplified form as 

- w"" - + + ^=^^r + (^-^ + 7.'»s) / + ^^^z = » (49) 
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where we take the vorticity n = 1 and ignore the contribution from A{r) given that this vanishes 
at the origin. Note that higher order terms in / coming from the potential are to be ignored in the 
hnearized analysis. We proceed to make a Frobenius ansatz for the behaviour close to the origin of 
the form 

oo 
m=0 

Upon substituting this ansatz in eq. ( |49| ) and looking at the lowest order in r one obtains the indicial 
equation of the system, hence we look at the equation at order r^~^ where we obtain 

(A-3)(A-1)2(A + 1) = 0. (51) 

Therefore we have three distinct roots A = 3,1,-1 with multiplicities 1,2,1 respectively. We 
proceed to determine the coefficients a^n by looking at higher orders in r. The equation at order 
j.A-3 ij]2piigs that = 0. The order r^~^ equation leads to 

which gives solutions for both roots A = 1 and A = 3 as 

oo oo 

/(i)=^E-^-(l)^"' /(3) = E -^-(3)r- (53) 

m=0 m=0 



where 02 and 62 are coefficients extracted from eq.(52) with the appropriate choice for A, and an 



ill-defined solution for A = — 1 which we will return to later. The solution /2 corresponds to the 



behaviour used in eq. ( 20 ) at n = 1 . Both these solutions and their derivatives are well behaved at 
the origin. The next order coefficients can be extracted from the order equation as 

^ _ -7^0 («^(1 + 3A + A^) - je^alm + 4A + A^)) 

/32(l + A)(3 + A)2(7+17A + 8A2 + A3) ^ ' 

with higher order J>n's for odd m vanishing. Being A = 1 a multiplicity 2 root, we know that the 
A derivative of this solution is also a solution of the equations of motion. In general if /(jj^) is a 



solution of the form eq.(50), then 



^ = In r/(,)+r ^ (55) 

m=0 

which means that an independent solution is of the form 

/(I) = ^ = ^ £ -^-(1)^" Inr + r f; ^^r™. (56) 

m=0 m=0 

This solution has a singular derivative at the origin. 

The solution of the linearized problem for A = — 1 takes the form 

^00 00 

m=0 m=0 

where as before the sum extends over even m and one finds that Bq and Cq are non- vanishing. This 
solution of the linearized problem diverges at r = and hence should not be taken into account for 
searching physically acceptable solutions. 

Acknowledgments : This work was supported by CONICET , ANPCYT , CIC, UBA and UNLP, 
Argentina. 



12 



References 

[1] H. B. Nielsen and P. Olesen, Nucl. Phys. B 61 (1973) 45. 

[2] V. L. Ginzburg, Pis'ma Zh Eksp. Fiz. 30 (1979) 345 (JETP Letters 30 (1979) 319). 

[3] A. I. Buzdin and M. L. Kulic, J. of Low Temp. Phys. 54 (1983) 203. 

[4] H. W. Diehl, Acta Phys. Slov. 52 (2002) 271. 

[5] P. Horava, Class. Quant. Grav. 28 (2011) 114012 |arXiv: 1101. 1081] [hep-th]] . 

[6] S. Kachru, X. Liu and M. Mulligan, Phys. Rev. D 78 (2008) 106005 |ar Xiv : 0808 . 17251 [hep-th] ] . 

[7] R. M. Honreich, M. Luban and S. Shtrikman, Phys. Rev. Lett. 35 (1975) 1678. 

E. Ardonne, P. Fendley and E. Fradkin, Ann. of Phys. (N.Y.) 310 (2004) 493 Phys. Rev. Lett. 
66 (1991) 276. 

[9] S. Deser, R. Jackiw, S. Templeton, Phys. Rev. Lett. 48 (1982) 975; Annals Phys. 140 (1982) 
372-411. 

[10] S. K. Paul and A. Khare, Phys. Lett. B 174 (1986) 420 [Erratum-ibid. 177B (1986) 453]. 

[11] H. J. de Vega and F. A. Schaposnik, Phys. Rev. Lett. 56 (1986) 2564; Phys. Rev. D 34 (1986) 
3206. 

[12] J. Hong, Y. Kim and P. Y. Pac, Phys. Rev. Lett. 64 (1990) 2230. 
[13] R. Jackiw and E. J. Weinberg, Phys. Rev. Lett. 64 (1990) 2234. 

[14] L. F. Cugliandolo, G. Lozano, M. V. Manias and F. A. Schaposnik, Mod. Phys. Lett. A 6 
(1991) 479. 

[15] C. -k. Lee, K. -M. Lee and E. J. Weinberg, Phys. Lett. B 243 (1990) 105. 
[16] G. Grinstein, Phys. Rev. B 23, 4615 (1981) . 



[17] P. Hofava, Phys. Lett. B 694,172 (2010) |arXiv:0811.2217| [hep-th]]. 



[18] M. Mulligan, C. Nayak and S. Kachru, Phys. Rev. B 82 (2010) 085102 | arXiv:1004.3570| [cond- 
mat.str-el]]. 

[19] A. Michelson, Phys. Rev. B 16, 577 and 585 (1977). 

[20] E. A. Coddington and N. Levinson, Theory of Ordinary Differential Equations, Tata McGraw- 
Hill Pub. Co., New Dehh, 1972. 



[21] M. Mulligan, C. Nayak and S. Kachru, Phys. Rev. B 82, 085102 (2010) |ar Xiv: 1004.3570 
[cond-mat.str-el]]. 



[22] L S. Landea, N. Grandi and G. A. Silva, arXiv:1206.0611 [cond-mat.str-el]. 



13 



